The reflection-antisymmetric counterpart of the 



Karman-Howarth dynamical equation 

Susan Kurien 

Center for Nonlinear Studies and Theoretical Division, 
Los Alamos National Laboratory, Los Alamos, New Mexico 

February 8, 2008 

Abstract 

We study the isotropic, helical component in homogeneous turbulence 
using statistical objects which have the correct symmetry and parity 
properties. Using these objects we derive an analogue of the Karman- 
Howarth equation, that arises due to lack of mirror-reflection symmetry 
in isotropic flows. The main equation we obtain is consistent with the 
results of O. Chkhetiani [JETP, 63, 768, (1996)] and V.S. L'vov et al. 



| http:/ /xxx. lanl.gov/abs/chao-dyn/9705016 , (1997)] but is derived using 



only velocity correlations, with no direct consideration of the vorticity 
or helicity. This alternative formulation offers an advantage to both ex- 
perimental and numerical measurements. We also postulate, under the 
assumption of self-similarity, the existence of a hierarchy of scaling ex- 
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ponents for helical velocity correlation functions of arbitrary order, anal- 
ogous to the Kolmogorov 1941 prediction for the scaling exponents of 
velocity structure function. 

1 Introduction 

In their 1938 paper on the statistical properties of homogeneous, isotropic, 
reflection-symmetric turbulence, T. von Karman and L. Howarth derived the 
equation for the dynamics of the two-point velocity correlation function 
This equation is of fundamental importance since it relates the mean rate of 
change of energy to the flux of energy across a given correlation length r in the 
flow. A form of this equation was used by A.N. Kolmogorov in 1941 |^ (K41) 
to derive one of the few exact results known for isotropic, homogeneous, and 
reflection-symmetric turbulence, the "4/5ths law" which relates the third-order 
longitudinal structure function to e, the mean rate of energy dissipation 

((ui(x + r)-u^(x))-^> = -ler (1) 

where U/, is the component of the velocity along the separation vector r. If 
the flow is not reflection-symmetric however, a new equation may be derived 
to complement the Karrhan-Howarth equation. Three recent works have de- 
rived equations for third-order statistics in isotropic helical flows by considering 
velocity- vorticity correlations |4| ||. In this paper, we show that the Karman- 
Howarth equation has a counterpart which arises due to parity-violation in 
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isotropic flows and which can be written solely in terms of two-point velocity 
correlations. We demonstrate the equivalence of our result with those of Q and 

i- 

We were motivated in this work by a series of investigations which proposed 
the use of the S0(3) decomposition of tensor quantities, the structure functions, 
defined by 

Saf}{r) = ((iia(x + r) - Ua{x)){up{x + r) -u^(x))) (2) 

in order to study the anisotropic contributions to their scaling. The decomposi- 
tion of the structure function into rotationally invariant, irreducible subgroups 
of the S0(3) symmetry group S-j^^ {r)+S^^^ {r) + ... allowed the separation of the 
isotropic (indexed by j = 0) from the anisotropic (indexed by j > 0) contribu- 
tions to the structure function. This procedure has allowed better quantification 
of the rate of decay of anisotropy of the small scales in turbulence |^ • These 
analyses considered homogeneous, isotropic and reflection symmetric flows. In 
the isotropic (j = 0) sector, the reflection symmetric structure function tensor 
has the form 

= C,ir)d^f3 + C2{rf-^ (3) 

Homogeneity and incompressibility provide a constraint between the scalar func- 
tions Ci(r) and C2(r). If the condition of reflection symmetry is dropped, there 
arises a further tensor contribution to the isotropic sector given by SafSj^- This 
contribution is interesting because it is isotropic (rotationally invariant), which 
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implies that it belongs in the j = sector, but is antisymmetric in (a, (3) and 
changes sign under mirror reflection of r. Since the second order structure func- 
tion is symmetric in its indices and does not change sign under inversion of 
r, it simply cannot to be used to observe this antisymmetric contribution. In 
fact, when the antisymmetric contribution is included in our decomposition, we 
are effectively using the isotropic irreducible representation of the 0(3) sym- 
metry group which includes operations that are not reflexion invariant under 
r — > — r. Said differently, the elements A of the orthogonal group 0(3) satisfy 
det{A) — ±1. The elements with determinant +1 form the S0(3) symmetry 
group of all (even-parity) rotations while those with determinant -1 are (odd- 
parity) reflections. The present work demonstrates how to access this isotropic, 
antisymmetric, odd-parity contribution using the tensor object with the ap- 
propriate parity and symmetry properties. The dynamics of such an object 
will provide the antisymmetric counterpart to the Karman-Howarth dynamical 
equation. 

In section 2, we present and discuss the second- and third-order velocity 
correlations and their symmetric and antisymmetric contributions. In section 
4 we derive the antisymmetric, odd-parity counterpart of the Karman-Howarth 
equation for the second-order correlation function and show its equivalence to 
previous results. In section 5 we postulate the existence of generalized helical 
higher-order velocity correlations and their scaling behavior under the assump- 
tion of self-similarity. Section 6 provides a summary and discussion. 
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2 The symmetry and parity properties of the 
two-point velocity correlation functions 

2.1 The second-order correlation tensor 

The two-point correlation tensor function of velocity fluctuations is defined by 

Ra/sir) = (wa(x)w/3(x + r)) (4) 

where r is the vector separation between two points, and subscripts a, [3 are com- 
ponents in a chosen Cartesian coordinate system. In homogeneous, isotropic, 
and not necessarily reflection-symmetric turbulence, the correlation function 
may be written as a sum of the dyadics flG, O] 



Such a tensor may be written as the sum of its symmetric (in a, /3) and 
antisymmetric components as 

Raf3(r) + Rpa(r) Raf3{r) - Rpair) 

Rapi^) = ^ H 2 

= i?f^(r) + i?^^(r) (6) 

The symmetric contribution R^^ (r) consists of the first two terms on the right 
side of Eq. (||) while the antisymmetric contribution i?^^(r) is the last term in 
Eq. (i. 



If the flow is statistically homogeneous, then the incompressibility constraint 
implies 

daRo^pir) ^ dpRa,p{r) ^ (7) 

where da{-) denotes the partial derivative with respect to r^. The incompress- 
ibility condition applies separately to the symmetric and antisymmetric compo- 
nents as daR^piy) = djjR^piv) = and daR^piv) = dpR^pir) = since the 
symmetric and antisymmetric components are of opposite parity. This is an in- 
teresting and useful property of these correlation functions in the isotropic sector 
and for homogeneous flows - decomposition into symmetric and antisymmetric 
components automatically separates the even- and odd-parity contributions. 
The symmetric part i?f^(r) with tensor basis as follows, 

i?f^(r)=Ai(r)<5„^-f A2(r)^ (8) 

has been analyzed extensively (see for example, |^) under the assumption of 
homogeneous, isotropic and mirror-symmetric turbulence. These three condi- 
tions imply the translational, rotational and reflectional invariance respectively 
of a given statistical quantity used to describe the flow. Note that the struc- 
ture function (Eq. (|^)) is twice the symmetrized correlation function R^^ plus 
twice the mean-square velocity fluctuation. The latter addition makes the struc- 
ture function galilean invariant and hence a suitable candidate for the study of 
universal statistics of the small scales. 

The form of the antisymmetric tensor in the j = sector of the 0(3) repre- 
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sentation is 



-^a/3(r) = (-Ua(x)M;3(x + r)) - (u/3(x)ua(x + r)) 



= H{r)e^p^— (9) 



Let us apply the incompressibility constraint to the antisymmetric tensor form: 



r H(r) 
dH{r) 

^ (10) 



In going from the second to the last lines of Eq. (|T^, we have used the fact 
that contracting an antisymmetric tensor with a symmetric one gives identically 
zero. We conclude that incompressibility does not provide any constraint on the 
scalar coefhcient H{r) of the antisymmetric tensor contribution. 

We can give an argument that the antisymmetrized correlation function is 
galilean invariant by definition. Suppose we are in a frame moving with velocity 
U, then 

Ripir) - (K(x) + U)(m0(x + r) + U)) - ((u^(x) + U)K(x + r) + U)) (11) 

It is seen that, because of homogeneity and the minus sign used to antisym- 
metrize, any dependence on U drops out. Therefore, we can hope that, as in 
the case of the structure functions, the object i?^^(r) will display the (universal) 



properties of the small scales. 

2.2 The third-order correlation tensor 

Our goal is to derive the dynamical equation for the second-order antisymmetric 
correlation R^^ (r) as a counterpart to the Karman-Howarth dynamical equation 
for the second-order symmetric correlation i?^^(r) (denoted in their paper of 
1938 by Rik{^))- Since such an expression will involve the two-point third- 
order correlation function, we will first review its properties. 



has the following properties in homogeneous turbulence. It is clearly symmetric 
in indices a, 7, with mixed symmetry in other combinations a, (3 and 7,/3 and 
in general of mixed parity. By "mixed" we mean that the symmetry and parity 
properties are indeterminate. In the isotropic tensor representation then, there 



(12) 



are four terms [|[ || 




+ S{r){ 



(13) 



In anticipation of separating the terms of opposite symmetry as was done in the 
case of the second-order correlation function, we write 
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where S^^ ^ is antisymmetric in a, f3 and has tensor contributions as follows 
. _ (ua(x)M-y(x)M;3(x + r)) - (u^(x)M^(x)uQ(x + r)) 

'^Q7,/3ll"j — ^ 

Si{r) - S2{r) ( rp . 

= ^ [S^f--^i^-fy) 

These are the terms which were excluded in the Karman-Howarth equation for 
reflection-symmetric flows. 



3 The antisymmetric component of the Karman- 
Howarth equation 

We now derive in a simple manner the dynamical equation for i?^^(r). As in 
Hinze's |^ equation 1.48, starting from the Navier-Stokes equation for homoge- 
neous turbulence we can write the equation for Rap 

d 1 

-Q^Rap — drjSa^^p + d^Sa,jj3 — — -( — l9a^p,/3 + dpKa^p) + 2l^dj^Raf3 (16) 

where Ka.p = {ua{'x.)p{'x.+r)) andp is the pressure. We write a similar equation 
for Rpa which we subtract from Eq. ( p^ and divide throughout by 2. 

d / Ral3 — Rpa \ „ / Sg^^p — Spy, a ^ / Sg^^p — Sp,ya \ 

dtX 2 J ~ 2 J+'^n 2 J 
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^ ( - ^{-daKp^p + dpKa,p) + ^{-dpKp^a + daKp^p)^ 

2vd^^[^ ^ ) (17) 



The pressure terms may be shown to vanish identically using homogeneity 
and incompressibility and assuming regularity as r ^ 0, as in the reflection- 
symmetric, isotropic case The homogeneity condition Sa-y.p{r) = 
'S'/3.7q(— r) adds a further constraint, giving 



d 

'Ol^ap - '^djS^^ p = 2vd.f^R^p. (18) 



This equation is the antisymmetric counter-part to the Karman-Howarth equa- 
tion for the second-order correlation functions. All the quantities in this equa- 
tion are relatively easily measured in experiments and numerical simulations 
since no velocity derivatives are involved in the correlation functions, only the 
velocities themselves. Substituting in Eq. (Hq) the tensor forms for the anti- 



symmetric correlation functions (Eqs. (|9|) and (15)) we arrive at the dynamical 
relation relating the scalars H{r) and S{r) 



This equation was derived by Chkhetiani using the dynamics of velocity- 
vorticity correlations. In the present derivation, we have arrived at the conclu- 
sion without the need to directly consider vorticity or helicity. We only used the 
0(3) tensor representation for the correlation function in homogeneous, isotropic 
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flows in which symmetry and parity properties are trivially separated. 
3.1 Derivation of KH-helical scaling law 

We apply the curl operator to the second-order antisymmetrized correlation 
fmiction Eq. and obtain the leading order behavior of H{r) ~ Ti.r/3 (see 
Eq. ^ and associated details in the Appendix) where the mean hclicity Ti. = 
(u • (j^>)/2. We now substitute this leading order dependence of H{r) back into 
the KH law, 

l(¥-'0-(^l-?)-«-KI^-7|-^)-« ™ 

Here, if we make the same assumption as that the main contribution to 
the time-derivative comes from the linear term with the next order terms not 
changing in the inertial range, and neglect the right-hand side in the limit as 
1/ ^ 0, 

Sir) ^ (21) 

where h is the mean helicity dissipation rate. This agrees with the scaling law 
derived in (There is a difference of a factor of 1/2 in the definition of mean 
helicity between Q and the present work.) The assumption made in deriving 
this law is that we have fully developed, freely decaying turbulence. These are 
the same assumptions made by Kolmogorov in deriving the 4/5ths law and the 
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energy spectrum. It is with this assumption that the foUowing holds 



^ iy{{dkv^)idku;^)) = h. (22) 



If a driving force is introduced, additional terms arise in the helicity balance 
equation [2^ (for example {f ■ lu)) which may not directly allow us to derive Eq. 
(^l|). It is however, not unreasonable to expect that Eq. (|2^) will hold for the 
steady-state, forced high-Reynolds number case. An argument similar to that 
which Frisch used to prove the 4/5ths law for the forced case, would have 
to be used. This aspect will not be covered in the present work. 

In the Appendix we show that an alternative form of Eq. ^ may be derived 
in the form of the following pair of equations 



UL • (UT X O = -^/ir^ (23) 
UT ■ {{ul X Wrp) + (ut X u'^)) = -^hr'^ 



where the velocity vector has been separated into its longitudinal (along the 
separation vector r) and transverse components as u = + uj-. The un- 
primed velocities denote their value at x and primed velocities denote their 
value at X -f r. The first line of Eq. ^ is equivalent to the so-called "2/15ths 
law" derived by L'vov et al B (see Appendix for more details). 
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3.2 The scaling behavior of higher-order correlation func- 
tions 

The antisymmetrized correlation functions may be thought of as newly defined 
structure functions appropriate for helical flows. In the second and third order 
(i?^^ and S^^ p respectively) we have shown that the antisymmetrized correla- 
tion functions are galilean invariant, so that sweeping effects are eliminated as 
in the case of the symmetric structure functions, and we may hope for universal 
properties for small r. For the third-order correlation, we have seen that the 
scaling in the inertial range is ~ (Eq. (pTf)). Let us now make the K41 
assumption of self-similarity such that S^^ ^ ^ (-^a/?)"^^^ obtain the inertial 
range behavior for i?^^ ~ r^^^. This corresponds to an inertial range scaling 
of the (helicity) cospectrum £^12(^3) ^ k^'^/'^. The denotes the wavenumber 
component in the direction mutually orthogonal to a = 1, and (3 = 2. This 
estimate for the scaling of the cospectrum coincides with the Lumley 1967 esti- 
mate |l^ for the (anisotropic sector j = 2) shear-stress (Reynolds) cospectrum 
Ei2{ki). However, the present dimensional estimate for the cospectrum -©12(^:3) 
is due to the reflection symmetry breaking, not due to the rotational symmetry 
breaking. 

If we construct nth-order antisymmetrized correlation functions with scaling 
exponents in the inertial range, the self-similarity argument dictates that 
^„ = This would be the helical counter part to the K41 scaling prediction 
for the structure functions which says that the nth-order structure functions 
have scaling exponents Cn f • It is not at all clear that self-similarity is a 
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reasonable assumption to make even in the case of low-order helical statistics 
[ p7| . This conjecture may only hold in the case of the maximal helical cascade 
in which there is no joint cascade of energy p8|. 



4 Conclusion 

The understanding of helicity dynamics in three-dimensional flows is still evolv- 
ing. It has been known for some time that helicity is conserved in the fluid 
equations in the inviscid limit [ p^ . The simultaneous existence of both helic- 
ity and energy cascades to the high-wavenumbers was first considered by A. 
Brissaud et al. pc| |. In that work, the scenario for a pure helicity or maxi- 
mally helical cascade was also proposed, in which energy cascade to the small 
scales is blocked, giving rise to an energy spectrum E{k) - k-'^/^. R.H. Kraich- 



nan showed |21 , based on physical considerations, that the scenario of joint 
energy and helicity cascades to the high-wavenumbers, with recovery of the 
Kolmogorov energy spectrum E{k) ^ k^^/'^ is more likely. This joint-cascade 
picture has subsequently been strengthened by observations in numerical sim- 
ulations from which it seems likely that the helicity injected at the large 
scales cascades downscale, more or less passively transported by the energy cas- 
cade. More recently, Ditlevsen and Giuliani show, both theoretically |Q and 
using shell-model calculations |Q, that at high- Reynolds numbers a joint cas- 
cade of energy and helicity must exist in some range of wavenumbers. They 
argue that for wavenumbers larger than this range the reflection-symmetry is 
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restored by the dominant helicity dissipation term. Q. Chen et al. |^ have 
shown by means of heUcal-wave decomposition of the velocity field, that the de- 
tailed transfer of energy (and helicity) between helical-wave modes of opposite 
parity is consistent with the existence of a joint cascade, with —5/3 scaling for 
both energy and helicity spectra. They also confirm their theoretical predic- 
tions using numerical simulations. Their analysis disagrees with p3| , over 
the precise range wavenumber over which these cascades exist, but nonetheless, 
both works agree that for high Reynolds numbers, a joint cascade of energy and 
helicity will coexist for some range of wavenumbers, with parity restoration at 
sufficiently large wavenumbers. The present analysis is also consistent with the 
joint cascade at high Reynolds numbers. The original Karman-Howarth result 
and the helical version derived here are not mutually exclusive. The former 
picks out the reflection-symmetric part of the flow, while the latter picks out 
the reflection- antisymmetric part. The two contributions are measured by dif- 
ferent quantities which allow for the separation of the parity and symmetry 
properties of the flow. It is not clear, atleast to this author, whether the present 
formulation predicts that energy and helicity cascades will coexist for all scales 
(consistent with ||2^) or for only a certain range of scales (consistent with |p3| ) 
at Reynolds numbers high enough. This might have to be left to empirical tests 
using experimental and direct numerical simulations data. Thus far, only the 
shell-model simulations of Ditlevsen and Giuliani and Biferale et al 
exist to guide our intuition as to the scaling ranges of the cascades. 

From the analysis of pa it appears that if one of the helical modes is blocked, 
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which can easily be done in simulations and shell- model calculations, but may 
not be possible in real flows, a pure helicity cascade will develop in the remaining 
mode which blocks the energy cascade down to small scales, and yields an energy 
spectrum E{k) ^ fc~^/'^. In this sense, the dimensional (self-similarity) argu- 
ment of section 3 for the scaling exponent of the cospectrum £^i2(fc3) '^^^ is 
consistent with the scenario of a pure helicity cascade. Without speculating on 
the feasibility of physically achieving such a purely helical cascade, we remark 
that the cospectrum £^12(^:3), of the two orthogonal components along the third 
orthogonal direction, is a fundamentally different object than the energy spec- 
trum and may well display entirely different functional behavior. The helicity 
cospectrum, which should be identically zero for homogeneous, isotropic, reflec- 
tion symmetric turbulence is a sensitive measure of reflection symmetry breaking 
[ p7| and the presence of helicity. The present work has shown that it arises from 
precisely that contribution to the second-order correlation which was excluded 
in the original isotropic, homogeneous, reflection-symmetric Karman-Howarth 
equation. 

A further new possibility suggested by this work is the construction of anti- 
symmetric higher-order (greater than 3) correlation functions. Assuming each 
pair of indices of an nth order, two-point velocity correlation function can be 
appropriately antisymmetrized as has been done here for the second and third 
order cases, we may have a new series of objects, which we will call helical 
structure functions. These, along with the usual structure functions familiar 
from studies of isotropic, reflection-symmetric flows, would form a complete set 
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of statistical objects with which to investigate statistical turbulence theories 
which are not necessarily confined to reflection-symmetric configurations. All 
the usual issues such as scaling exponent values, intermittency and anomaly 
could be studied for the helical structure functions. This work provides a di- 
mensional argument for what their scaling exponents could be. It is of interest 
to see how these behave relative to our predictions and to the anomalous scaling 
known for the usual structure functions. 

In conclusion, the present approach taken to derive the antisymmetric, or 
helical Karman-Howarth equation is not inconsistent with other recent work 
[| |, |. However, our derivation directly studies the dynamics of precisely 
those components of the second-order correlation functions that were omitted 
in the Karman-Howarth equation because of assumed non-helicity of the flow. 
The information about helicity of the flow is then obtained from velocity corre- 
lations instead of velocity-vorticity correlations or other correlations involving 
velocity derivatives. In a high-Reynolds number experimental flow, measuring 
the second-order antisymmetric velocity correlation function R^2{^) ^ coor- 
dinate system chosen such that r = rk, would give information about the mean 
helicity in the flow, while measurement of Sf^2 (^) would give information about 
the helicity flux. Such objects are ideal candidates for detecting parity violation 
in flows without having to resort to direct measurement of helicity. We intend 
to present the related analysis of numerical and experimental data in a future 
work. 
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A Appendix 

The antisymmetric tensors i?^^ and S^^ ^^ are directly related to helicity dis- 
sipation and fluxes as we will now demonstrate. Let us consider a particular 
geometry in which the separation vector r is along the 2(3)-axis. Then the only 
non-zero components are R12 = —R21 = (^ii(x)u2(x + r) — 'U2(x)ui(x + r)). 
This particular object is the correlation of the two components of the velocity 
orthogonal to the vector r. In the usual convention, it is the correlation of the 
V and w components relative to the separation vector. The only non-zero con- 
tribution to the isotropic, antisymmetric velocity correlation tensor comes from 
the velocity components orthogonal to the separation vector. 
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u,(x) 




Figure 1: Caricature of the type of correlation functions i?^^ which are non- 
zero in flows that are not reflection-symmetric. The curved arrow indicates the 
"handedness" of the fimction. 

We contract the tensor R^^ with the antisymmetric tensor eap^y p^ . 

eaf3j{Ua{-X.)up{:X. + r)) + ega-y(M^(x)Ma(x + r)) 

2 

(u(x) X u(x + r))^ 

The (pseudo)scalar function H{r) is the mean cross product of the velocities 
at two points separated by the vector length scale r. To choose a particular 
coordinate system, if the separation vector r lies along the z-axis, then H{r) is 
the z-component of the the cross-product. It vanishes as \r\ 0. The physical 
picture is depicted in the cartoon of Fig. 0. This result is to be compared with 
the corresponding result for the symmetric contribution R^^ contracted with 
6ap. In that case, what is obtained is {ua{x)ua{x + r)) = (u(x) • u(x + r)), 
the mean scalar (dot) product of the velocities at two points separated by the 
scale r; as r — > 0, we recover the non-zero mean energy ~ (u^). 



f-al3iH[r)eal3u — 

r 

^ 2H{r)^. (24) 
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The function Ri2{f3) may be thought of as a measure of momentum transfer 
between two orthogonal components of velocity along the direction perpendic- 
ular to both of them. If we take the curl of R'^/^, we have 



(u(x)..;(x + r)) . 2^+4^. (25) 

or r 



Taking the limit as r 0, the left hand side reduces to (u • a;) = 2'H where Ti, 
is the mean helicity of the flow, and we can solve for what must be the leading 
order behavior of H{r). 



H{r) = \nr + ... (26) 



The scalar coefficient H{r) of the antisymmetric tensor i?^^(r) is, in the leading 
order, a direct measure of the mean helicity in the flow. If we consider the 
particular coordinate system of Fig. |^ we see that Ri2{'r) = H{r) which is 
a leading order measure of the mean helicity of the flow. We note again the 
advantage of this formulation which allows one to measure mean helicity using 
only velocity correlations, without having to measure any local gradients. 

We perform a similar analysis for the third-order object with the contraction 

(z.,(:r)(u(x)xu(x + r))^) = 5(r) - <5,^) (27) 
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If we now proceed to write, as in Q, the velocity vector as the sum of its 
longitudinal (along r) and transverse components such that u(x) = Ui(x) + 
ut(x), we have 



(K+ut)^(K + ut) X K+u^r))^) = 5(r)(3^-(5^^) (28) 

where the un-primed velocities denote measurement at x and the primed veloc- 
ities denote measurement at x + r. It is clear that both the left and right side 
vanish for 7 = /i. But we would like to examine the detailed balance in terms 
of the longitudinal and transverse components on the left hand side. 

M^^ = (^{UL + Ut)^ ( {ul X U^)^ + {ul X Uy)^ + (ut X u'^)^ + (ut X u^)^ 



-S{r)[3^^S,,). (29) 

Since this must be true for any choice r, we can, without loss of generality, 
choose r = ri so that it lies along the x-axis. The matrix of Eq. (^) is then 
diagonal and traceless, and we see, using Eq. Gl^, that 



Mil = Ul • (ut X u^) = 2S{r) = ^hr^ (30) 
M22 = M33 = ut ■ {(ul + ut) X (u'l+u't)) ^ -Sir) ^ -^hr^ 



Eq. (gO|) is a form of the so-called "2/15ths law" derived by L'vov et al fj]. The 
exact result of that work is obtained by computing (u^ — u^) • (u^ x u^)) which 
is equal, by homogeneity to 2ul ■ (uy x u'rp) = 4S'(r) — j^hr'^. 
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